Searching for Smoking Gun Signatures of Decaying Dark Matter 



Joshua T. Ruderman 1 and Tomer Volansky 2 

'Department of Physics, Princeton University, Princeton, NJ 08544 
2 School of Natural Sciences, Institute for Advanced Study, Princeton, NJ 08540 

Clear methods to differentiate between decaying and annihilating dark matter (DM) scenarios are 
still by and large unavailable. In this note, we study the potential astrophysical signatures of a new 
class of hidden sector decaying DM models, which can address the recent cosmic ray measurements. 
Such models may produce primary photons and/or neutrinos at large rates, correlated with the 
leptonic production. The photon and neutrino spectra will then contain sharp features at the TeV 
scale. We demonstrate the discovery potential for upcoming and future measurements by FERMI, 
HESS, AGIS and IceCube/DeepCore. We show that these models may be discovered in the near 
future. Specifically, measurements of diffuse gamma rays by FERMI can detect the start of a hard 
photon feature. We argue that these hard spectra can be produced by decaying dark matter and be 
consistent with current constraints, but are difficult to reconcile with models of annihilating DM. 
Consequently the measurement of a hard spectral feature, in correlation with the current cosmic ray 
measurements, will strongly favor decaying DM models. Finally we comment on the preliminary 
results from the Inner Galaxy presented by the FERMI collaboration. 



Introduction. Recent astrophysical measure- 
ments [1] may constitute indirect detection of dark mat- 
ter (DM), the validation of which still awaits unam- 
biguous evidence. Such evidence may take the form 
of associated signals in photon and neutrino channels, 
which accompany the leptonic excess. If produced from 
DM decays or annihilations, the energetic leptons imply 
that photons are produced from Bremsstrahlung radi- 
ation (FSR) and Inverse Compton Scattering (ICS) off 
starlight, dust, and the CMB. There can also be neutri- 
nos or photons from lepton or pion decays. Even if these 
associated signals are discovered, the nature of dark mat- 
ter will remain, in part, obscured. In particular, it will 
be hard to determine whether an excess arises from de- 
cays or annihilations of dark matter. The reason for this 
is that the decay and annihilation fluxes scale differently 
with the density. Consequently, the flux differs between 
the two scenarios mainly at the Galactic Center (GC), 
where the density is highest, but also where the back- 
ground dominates. 

A more promising situation may arise if sharp spec- 
tral features are measured. Such features have not been 
conclusively observed in the electron spectra, but may 
appear in photon or neutrino measurements. For the re- 
mainder of this letter, we loosely define the hard photons 
or neutrinos that contribute to a sharp spectral feature, 
as 'primary'. If discovered, this feature would be smok- 
ing gun evidence for DM. Here we point out that such a 
signal, if found to be correlated with the leptonic excess, 
also has the potential to differentiate between the decay- 
ing and annihilating scenarios. Indeed, annihilating DM 
models that address the leptonic excess and have a sizable 
branching fraction into primary photons or neutrinos, are 
already excluded assuming the DM profiles supported by 
current N-body simulations [2 . Conversely, in the de- 
caying case, the photon or neutrino signals lie just below 
the sensitivity reach of current experiments. As we show, 



these may be measured in the near future thereby not 
only confirming the presence of DM, but also supporting 
the decaying DM scenario. 

The interpretation of the cosmic ray data requires the 
DM to have mass ~ few TeV. Many models have been 
put forth to explain the data (see [5] and refs. therein) 
nearly all of which bifurcate into annihilating and decay- 
ing DM scenarios. In order to explain the measurements, 
models of annihilating DM require a surprisingly large 
cross-section, three orders of magnitude larger than that 
of a thermal weakly interacting massive particle (WIMP) . 
Such an enhancement can, in principle, be achieved at 
low velocities, for example, by the Sommerfeld effect [4], 
without altering the compelling features of the WIMP 
scenario. In practice, however, it may be difficult to ob- 
tain a large enough enhancement and even if achieved, 
annihilating DM models are in tension with various 7- 
ray and neutrino bounds [3j [5] , and constraints from the 
recombination epoch [6]. 

The alternative of decaying DM solely replaces the 
need for the Sommerfeld enhancement. If dimension six 
operators are generated at or around the GUT scale [7J, 
the decay of the DM particle naturally explains the ob- 
served signals, again without altering the predictions of 
the thermal WIMP scenario. In this case, the various 
constraints from 7-rays, neutrinos and the like are sig- 
nificantly weaker and such signals may be probed by a 
variety of experiments [H] . The drawback of many decay- 
ing DM models lies in the difficulty to naturally explain 
the lack of hadronic activity, without introducing signif- 
icant fine-tuning or complicated and ad hoc structure. 
We present a new class of models where weak-scale DM 
decays into a hidden 'dark' sector which is broken at the 
GeV scale and which communicates with the SM through 
gauge kinetic mixing [9]. The lack of hadronic activity 
is explained through kinematics, much as in the model 
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Below we study one representative example of this class 
of models and show that DM may naturally decay into 
primary photons and/or neutrinos with large branching 
fractions. We study the current and prospective mea- 
surements of the spectra at FERMI, HESS, AGIS, Su- 
perK and IceCube/DeepCore. Our best-fit of the model 
to the cosmic ray measurements is consistent with cur- 
rent constraints and within the sensitivity reach of the 
above experiments. We comment on the recent prelimi- 
nary FERMI results in the diffuse photon spectrum which 
are consistent with the predictions of our model. 

A Model. Let us briefly describe a natural model 
of decaying DM which predicts correlated high-energy 
gamma rays and neutrinos. To avoid overproduction of 
antiprotons, the DM must predominantly decay into light 
leptons. This can be achieved cither through additional 
symmetries that forbid couplings to hadrons or through 
kinematical constraints. The former is more difficult to 
accomplish, especially for decaying DM where the sym- 
metries must be present at the GUT scale, and typically 
requires some amount of fine-tuning. We concentrate on 
the case where weak-scale DM decays into a light state 
that subsequently decays into leptons, with antiprotons 
suppressed by kinematics |llj . The DM itself may or 
may not be charged under the SM gauge group. Here we 
take DM to be SM-charged, which results in a primary 
neutrino signal. Below we present a brief description of 
the model and refer the reader to [9] for a more compre- 
hensive analysis of this and similar decaying DM models. 

We work within the context of supersymmetric GUTs. 
The DM states, \ + Xi are taken to be two chiral super- 
fields with charges (5, 0) + (5, 0) under SU(5)s M x U(l) d 
where SU(5) S m => SU(3) C x SU(2) W x U(l) Y is the 
GUT gauge group. The predicted signatures to be dis- 
cussed below do not depend on the choice of GUT or dark 
gauge group. The U(l)d mixes with the SM through ki- 
netic mixing, 

- e -Jd 2 ew d w Y . (i) 

Here WdiWy) is the dark(hypercharge) field strength. 
This term can be generated by heavy bifundamental 
fields and e is naturally of order e ~ 10 -3 — 1CP 4 
[TOl fT2"] . The above gauge mixing induces an effective 
Fayet-Iliopoulos term proportional to e(Dy), naturally 
triggering the spontaneous breaking of the U(l)d at the 
GeV scale [T3] . If this is the dominant term for the break- 
ing, the dark sector is, to a good approximation, super- 
symmetric. Conversely, the U(l) d breaking can be trig- 
gered by SUSY-breaking effects communicated through 
TeV scale fields charged under the dark sector [ID] . In 
this case the dark spectrum is not supersymmetric and 
the dark gaugino can be heavier or lighter than the dark 
gauge boson. 



A dimension six operator of the form, 

i^/^Wj (2) 

induces the DM decay. Here a d = g d /^ is the dark 
gauge coupling and 5/ denotes a SM multiplet. Such 
an operator is generated by integrating out GUT-scale 
fields, 1,7 H|, 

W= (M GVT + X)YY + M GVT XX + X X 5f- (3) 

where X and X are singlets and Y + Y have charges 
(1,1) + (1,-1) respectively. There are other possible 
dimension six operators that can decay DM into the dark 
sector, and the choice of decay operator does not alter the 
qualitative conclusions of this work. The decay operator 
of Eq. ([3]) induces DM decay with a lifetime of order, 

T ~ 10 26 sec ( -^-Y 2 ( H^LY 5 ( V (A) 

which is the right time scale to explain the astrophysical 
anomalies. We stress, however, that the DM lifetime is a 
free parameter of the theory. 

As a consequence of Eq. ^ , DM decays into a neutrino 
or sneutrino and two dark gauge bosons or gauginos. The 
production of primary neutrinos is a generic consequence 
of DM being electrically neutral while carrying hyper- 
charge [13] . The other possibility of decaying into the 
neutral Higgs is disfavored due to the antiproton bound 
from PAMELA. The subsequent decay of the sneutrino 
depends on the MSSM spectrum. Here we assume that 
the sneutrino is the (SM) NLSP and can decay into a 
gravitino and a neutrino. If the sneutrino is heavier, its 
decay may produce antiprotons. The resulting antipro- 
ton spectrum, however, is expected to be rather soft due 
to the cascade decay. The production of the dark gauge 
bosons and gauginos is followed by their decay through 
the kinetic mixing. The gauge bosons decay into light SM 
leptons, 7(j — > l + l~ . To avoid the tension with the an- 
tiproton and gamma-ray bounds from the GC, the gauge 
boson mass, m 7d , is required to be sufficiently small, be- 
low ~ GeV. The gauginos, on the other hand, decay to 
either dark gauge bosons or photons, as we describe in 
the next section. 

There are several constraints on the model described 
above. If DM couples elastically to the Z, the model is 
ruled out by direct detection [15] by 2-3 orders of magni- 
tude |16j . This constraint can be evaded by introducing 
a dark matter splitting, <5mDM ^100 keV, such that DM 
couples inelastically to the Z. With a splitting, there is 
a strong constraint from neutrino telescopes on inelastic 
capture in the Sun, followed by annihilation into ZZ or 
W + W~ [17]. This constraint can be evaded by taking a 
slightly larger splitting of i5todm ^ 500 keV. We evade 
these constraints by coupling DM to the SM Higgs: 

Wspiit = m N N 2 + m DM X2X2 + \2H d N (5) 
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where N is a singlet with weak-scale mass Mjv, and \i 
denotes the doublet component of %, where the DM re- 
sides. DM stability requires Mjy > A/dm, and integrating 
out N generates the required DM splitting. 

If the dominant DM annihilation channel is to SM 
gauge bosons, its mass is constrained to be of order 
1.1 TeV, in order to obtain the correct thermal abun- 
dance [16 . Such a low mass is inconsistent with the 
FERMI and HESS measurements. Interestingly, the 
same operator that generates the splittings to evade the 
direct detection bounds, Eq. (p)]), also opens a new an- 
nihilation channels into Higgses which can easily dom- 
inate the DM annihilation cross-section. Consequently, 
the mass of the DM is a free parameter. 

Another important constraint on this model is re- 
lated to the lifetime of the triplet partner of the DM, 
\3- There are strong constraints on colored particles 
with T3 > 10 17 sec because they form exotic atoms [15] . 
These constraints are evaded if the triplet partner decays 
through a dimension-5 operator, which we take to be in 
the Kahler potential, 

where s is a singlet. This operator decays \3 m T 3 ~ 
1 sec, but keeps the DM stable as long as todm < fn a < 
m X3 . That to X3 is heavier than todm is a generic conse- 
quence of RG evolution. 

We note that there are a number of operators that must 
not be present, such as dimcnsion-5 operators which in- 
duce prompt DM decays, and tree-level Yukawa couplings 
between DM and SM matter. All such operators are 
easily forbidden by symmetries at the GUT scale. Cos- 
mology also places nontrivial and interesting constraints 
on the spectrum and lifetimes of light fields in the dark 
sector. For a detailed discussion of both GUT-scale sym- 
metries and cosmological constraints on this model, we 
refer the reader to [3] . See also for discussions of the 
cosmology of light hidden sectors. 

Gamma-ray Signatures. The measurement of a 
photon line or sharp spectral feature is considered to be a 
smoking gun signature for DM. Typically, such features 
are model dependent and are associated with the high 
end of the spectrum, set by the DM mass. Interestingly, 
within the context of the observed electronic activity, a 
measurement of this kind would point towards the de- 
caying rather than annihilating DM scenario. Indeed, as 
was shown in [3] [5], annihilating DM models that provide 
explanations of the PAMELA and FERMI anomalies and 
which have sizeable branching fractions into primary pho- 
tons or 7r°'s are excluded. The exclusion holds under the 
assumption of DM profiles suggested by state-of-the-art 
N-body simulations. To ameliorate it, one would have to 
resort to currently unsubstantiated profiles (see however 
|20j ) which predict a significantly shallower or flatter den- 
sity, p, at the GC. A key observation is then: if measure- 



ments of photons that are correlated with the PAMELA 
and FERMI results indicate the production of primary 
photons or photons from ir° 's, then decaying DM models 
are strongly favored. 

Models of DM in which the DM decays or annihilates 
into a light hidden sector, predict gamma rays from sev- 
eral channels. Three irreducible sources of photons are 
FSR, ICS and synchrotron radiation. Concentrating on 
the high energy spectrum, the first two dominate, where 
FSR dominates at the high end of the spectrum while the 
ICS dominates at lower energies. Two additional sources 
of photons are those produced from ir° — > 77 decays and 
primary photons. In this letter we concentrate on the 
latter. Whether primary photons are produced depends 
strongly on the dark spectrum. Assuming a light grav- 
itino, if the dark gauginos are heavier than the dark gauge 
bosons, fn ld > m ld , the decay jd — * 7d G dominates and 
no primary photons are produced. If, on the other hand, 
the dark gaugino is degenerate with or lighter than the 
dark gauge boson, it decays through the kinetic mixing 
to a SM photon and a gravitino: 

l d ^l + G. (7) 

Measurement of the primary photons produced through 
this decay is thus an indication of a light dark gaugino (or 
more generally a light fermion) in the dark sector. Con- 
versely, the absence of such photons indicates a heavy 
fermionic spectrum, and in particular implies significant 
SUSY breaking in the dark sector. We stress that as 
opposed to the primary neutrino signal to be discussed 
in the next section, and which follows from the specific 
operator of Eq. ([2]), the primary photon signal is rather 
generic for models where dark matter decays or annihi- 
lates into a light hidden sector, and the signal probes the 
lower end of the dark sector spectrum. 

In the model discussed above, the \ supermultiplet is 
split, and we assume that the fermion plays the role of the 
DM. We use 3-body phase space spectra with flat matrix 
elements, which captures the model-independent part of 
the spectrum. The resulting primary photons arise from 
the two available decay channels: 

X^vjdld, X~*vidld- (8) 

and the photon spectrum depends on the mass of the 
sneutrino, which we keep as a free parameter. We begin 
by fitting our model to the PAMELA, FERMI, and HESS 
electron and positron cosmic rays spectra, as shown in 
Fig. [TJ;. We obtain the best-fit values: m x — 3.3 TeV, 
m 9 = 370 GeV, m ld = 400 MeV, and r x = 4 x 10 26 
sec. Consequently, the dark photon decays into electrons 
and muons with branching fractions (0.75,0.25) respec- 
tively. We use this best fit model for the gamma ray and 
neutrino analyses that follow. 

We compute contributions to the photon flux from the 
primary photons described above, ICS, and FSR. No 
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FIG. 1: The electron fit and photon flux predictions for our best-fit model, using the Einasto DM profile with a = 0.17. The 
right plot shows the fitted curve to the e + + e~ flux measured by FERMI and HESS [I]. The middle curve shows the predicted 
spectrum from the Galactic Ridge. The red data points represent the corresponding HESS measurement with 2a error 
bars. The dashed blue line is the primary photon component predicted by our model, the FSR contribution is too small to be 
plotted, and the green line is the total flux, including a background modeled by a power-law and fitted to the HESS data. The 
signal is consistent with the current HESS data and can be constrained by further statistics. The black data points illustrate 
the projected reach for the AGIS experiment [30] with 2a error bars. As can be seen, AGIS will decisively test this scenario. 
The left plot shows the diffuse photon prediction in the 0° < I < 360°, 10° < \b\ < 20° region. The blue points correspond to 
the recent FERMI measurement [24]. The dashed lines show the primary (blue), ICS (black), and FSR (red) contributions to 
the total flux, shown in green. As the FERMI data extends to higher energies it may be able to measure the turnover between 
the ICS and primary photon components of the signal. 



7r u 's are produced due to the lightness of 7^. To com- 
pute the ICS we propagate the electrons following the 
procedure described in [3 . Throughout the analysis we 
assume the MED propagation model [21] and Einasto 
DM profile [22], 



p( r ) = P& Ex P 



{{rlr s ) a 1) 



(9) 



with r s = 20 kpc, p G = 0.3 GeV/cm 3 and a = 0.17. 
Since we're studying decaying DM, the dependence on 
the choice of profile is significantly weaker than in the 
annihilating case. Our gamma ray results are relatively 
insensitive to varying a in the range found in simulations, 
0.12 — 0.2, and to interchanging Einasto with NFW [2"3"] . 
The flux docs not change by more than a factor of 3 
around the GC nor does it vary by more than 10% for 
the diffuse gamma signal in the regions studied below. 

We consider signals from both the center of the Galaxy 
and diffuse gammas. A measurement of the latter in the 
region 0° < I < 360°, 10° < |6| < 20° was recently re- 
leased by the FERMI collaboration [24. The data and 
the predictions of the model are shown in Fig. [TJi. The 
blue dashed line corresponds to the primary photon spec- 
trum resulting from the DM decays, the red dashed line 
corresponds to the contribution from FSR and the black 
corresponds to that from ICS. As can be seen, the pri- 
mary photons dominate the high end of the spectrum, 
starting at ~ 100 GeV. Since FERMI can measure the 
spectrum up to 300 GeV, it may only probe the tail of 



the primary photon contribution. Although it will be dif- 
ficult for FERMI to disentangle the primary contribution 
from the ICS, it is possible that it will have enough sen- 
sitivity to resolve the turnover in the spectrum. Finally 
we note that a gauge boson with a mass larger than the 
400 MeV presented here, can produce a more significant 
FSR signal. For complimentary studies of such signals 
see [2"5] . 

A remark is in order. As we mention before, the mea- 
surement of the hard spectrum is a strong indication of 
decaying DM. The normalization of the spectrum is such 
that if a signal resulting from annihilation is seen in the 
diffused gamma region, away from the GC, it must be 
over-produced at the center itself. This is the case, for 
instance, when DM annihilates into two r's which conse- 
quently decays into 7r°'s (see Fig. 8 of [5]). This state- 
ment is in sharp contrast to the decaying DM scenario 
since the photon flux depends on p 2 in the former case 
and on p in the latter. Nevertheless, it is conceivable that 
the DM density will turn out to be shallower in which 
case the GC signal can be hiding below the background. 
A correlated measurement in the two regions would then 
imply the existence of decaying DM or some coincidence 
regarding the DM density at the GC. 

Measurements from the GC and Galactic Ridge (GR) 
have been taken by the HESS collaboration and already 
place strong constraints on models of DM ^\ . The GR 
is more sensitive since is measures a smaller flux in a 
larger region. Furthermore, it is less sensitive to the DM 
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FIG. 2: Neutrino constraints and discovery reach. The right plot shows the 95% C.L. muon flux limit from SuperK |31) 
as a function of the opening angle around the Galactic Center. The green curve corresponds to our best-fit decaying DM 
model, while the blue curve illustrates the expected flux for an annihilating model with a large branching fraction of primary 
neutrinos. Our decaying model is consistent with the constraints, while the annihilating model is excluded. The left plot shows 
the neutrino flux as a function of the energy. The green curve corresponds to the predicted total signal and the dashed lines 
show the primary (blue) and soft (red) neutrino components. The data correspond to the AMANDA measurement 37 . The 
black lines indicate the expected 1- (dashed) and 3- (solid) year 5a discovery reach in IceCube/DeepCore [34]. We find that 
there will be sufficient sensitivity to detect the flux in the last two bins of the predicted signal. 



profile since the GC itself (a 0.1° cone around the center) 
is subtracted. We show this measurement and the model 
prediction in Fig. [TJd, where the error bars correspond to 
2a. The red dashed line is the signal and the green line 
corresponds to the total signal plus background, where 
the latter is found by fitting a power-law to the data. Our 
predicted signal is below the current HESS resolution. 
There is also a preliminary FERMI measurement from 
the GC [27) and we find that the predicted signal of our 
model is well below the data. 

In addition to upcoming HESS results with improved 
sensitivity, there are two future Cherenkov telescope ar- 
rays under consideration: AGIS [28 and CTA [25]. Both 
arrays are expected to be significantly larger than HESS, 
with an order 5-10 improvement in sensitivity at the ener- 
gies of interest. To demonstrate the capabilities of these 
experiments, we concentrate on AGIS and project a pos- 
sible future measurement using an effective area of 1 km 2 , 
energy resolution of 15% and 200 hours of measurement. 
We show the projected measurement on Fig.[T]D, in black. 
The signal is found to be well above the expected sensitiv- 
ity |30j . To conclude, we remark that the future location 
of AGIS will influence its ability to take measurements 
from the GC. It is therefore preferable to locate it in the 
southern hemisphere, such as the proposed site near the 
Felix Aguilar Observatory in Argentina. 

Neutrino Signatures. The second smoking gun sig- 
nature of DM is a neutrino line or a sharp spectral fea- 
ture. While SuperKamiokande already places important 
constraints on DM models, the upcoming experiments 
IceCube, Antares and KM3NeT will have the potential 



to resolve features in the neutrino spectrum at the TeV 
scale and therefore discover DM. Here, by a sharp fea- 
ture, we mean a feature which appears in one or a few 
measured bins. In our model (or in corresponding annihi- 
lating DM models) , we therefore distinguish between pri- 
mary neutrinos (which are accompanied by a sharp fea- 
ture), produced directly from the decay or annihilation, 
and soft neutrinos that result from lepton and pion de- 
cays. Within annihilating DM scenarios that explain the 
lepton anomalies, primary neutrinos with a large branch- 
ing fraction are already excluded. On the other hand, as 
we show below, the corresponding decaying scenarios are 
right below the sensitivity of current experiments and can 
easily be discovered by future ones. 

As a first step, we must make sure that our model 
is not already excluded by the current neutrino bounds 
from SuperK. In Fig. |2]d, we plot the neutrino flux of 
our best fit model as a function of opening angle around 
the GC, integrated over energies detectable at SuperK. 
We find that it is consistent with the bounds of [3"T] . 
We note that annihilating DM models with order one 
branching fractions to primary neutrinos are excluded by 
this bound by 0(5 — 10). Therefore, the detection of a 
primary neutrino signal can favor decaying DM. This is 
illustrated in the figure where we plot the same spectrum 
for the case where the dependence on the rate depends 
on the density squared, as in annihilating DM. 

SuperK can also specifically identify primary neutrinos 
because at energies E v > 1 TeV, there is a high proba- 
bility of inducing electromagnetic showers. For our best 
fit model we calculate the muon flux from neutrinos in 
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a cone of 5° around the GC. The number of upward 
showering events is given by the expression |32j : 

Shower = J dE v ^J{E v )e{E v ). (10) 

Here <p v is the neutrino flux, / (E v ) is the conversion 
probability of neutrinos into muons of detectable energy, 
and e (E u ) is the showering probability extracted from 
the analysis of [32 . The result for our best fit model is 

Shower = 1-8 xlO" 16 cm" 2 s" 1 . (11) 

This is smaller than the expected background flux of 
atmospheric upward showering neutrinos, N£ = 3.4 x 
10~ 16 cm~ 2 s _1 . SuperK has done an analysis on up- 
ward showering neutrinos from 1646 days of data, and 
they find no events in the above cone [33] • Our best fit 
model is therefore not excluded, but interestingly close 
to the current sensitivity. 

The hard neutrino signal predicted by our model will 
be tested by upcoming neutrino experiments, such as 
IceCube and Antares. Here we focus on IceCube with 
DeepCore, which will have good sensitivity for TeV-scale 
neutrinos within the DeepCore fiducial volume |34j . The 
surrounding IceCube strings will be used to veto down- 
going muons, such that the main background is due to 
atmospheric neutrinos. The signal is then expected to be- 
come comparable to the background at high energies, be- 
cause the neutrino-nucleon scattering cross-section rises 
as <j u n ~ E 2 while the atmospheric neutrino background 
drops rapidly as E~ 3 . We estimate the future sensitiv- 
ity of IceCube with DeepCore by following an approach 
along the lines of |35j . The number of neutrino events to 
be observed at IceCube, as a function of neutrino energy 
is given by the expression, 

% = ^ P^ c N A a vN {E) V{E) Us , (12) 

where p icc = 0.9 g/cm 3 , N A = 6.022 x 10 23 g~\ i obs 
is the run time, <j v n is the neutrino-nucleon scattering 
cross-section, and V(E) ss 0.04 km 3 is the effective vol- 
ume of DeepCore for neutrino showers, roughly estimated 
from 34J . 

We compare our predicted signal at IceCube and Deep- 
Core with the atmospheric neutrino background [55] . In 
order to estimate the energy resolution of IceCube, we 
bin the background in bins of size log(-E max / E m i n ) = 0.3. 
In Fig.|2ji, we plot the AMANDA data [37] together with 
the flux sensitivity for 5a discovery in each bin, collect- 
ing events from the entire 2ir sky for 1 and 3 years of run 
time. We also plot the signal flux for our best-fit model. 
A 5cr discovery is found to be possible at the highest 
bin of the predicted spectrum, above E v > 1 TeV, af- 
ter 1 year, and in 2 bins after 3 years. This is in sharp 
contrast to the soft neutrinos which populate lower en- 
ergy bins and therefore cannot be discovered with this 



method. We see that IceCube with DeepCore has the 
potential to decisively test decaying DM with a primary 
neutrino component. It would be interesting to also es- 
timate the reach of Antares for this scenario. 

Note Added. While this work was in completion, 
new but very preliminary data from the Inner Galaxy re- 
gion, 0.25° < |6| < 4.75°, 0.25° < |/| < 29.75°, was pre- 
sented by the FERMI collaboration [38 . Interestingly, 
the results indicate a deviation from the naive expected 
background at high energy. Being preliminary, these data 
may well suffer from significant cosmic-ray contamina- 
tion, especially at the high end of the spectrum. Nonethe- 
less, we cannot resist to compare our predictions with the 
measurement, as we do in Fig. [3] In the figure, we plot 
two sets of curves for two different fits. The first (thick 
lower green line) is the best-fit model described above 
and the second (thin upper green line) corresponds to a 
lighter DM mass, resulting in a somewhat worse fit to the 
HESS electronic measurement. The ICS, hard and FSR 
contributions to the spectrum of the best-fit model are 
shown in dashed lines. A careful background analysis is 
warranted. Here we simply assume the dominance of 7r° 
decay at energies above ~ 1 GeV which can be approxi- 
mated by a power-law spectrum. We therefore model the 
background by fitting a power-law to the data between 3 
and 16 GeV. 

With the above cautious remarks in mind, we note the 
following: 

• If confirmed, the excess may indirectly imply the 
existence of DM [55] . 

• For our best fit model, with decay rate normalized 
to explain the PAMELA, FERMI, and HESS lep- 
tonic data, there is insufficient ICS to account for 
the excess at high energies > 100 GeV. The FSR 
that we find is also insufficient to account for the 
excess. 

• For our best fit model, the primary photon com- 
ponent is consistent with the preliminary data 
at high-energies. Still, our model's prediction is 
slightly below the data. We find two sources of ten- 
sion: (i) fitting to the HESS electron measurement 
results in a high DM mass, setting the hardness of 
the photon spectrum, and (ii) the decay rate can- 
not be increased without introducing tension with 
the HESS GC/GR measurements. We expect the 
overall fit to improve should the data points move 
down following a complete analysis. 

• Future improvement of the statistics may be capa- 
ble of indicating a turnover in the spectrum. 

• Ref. [40] presents an interesting analysis of the data 
and its implications to models of annihilating DM. 
We comment that the authors utilize a different fit- 
ting procedure from the one presented here. While 
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Diffuse Photons 




1(T 2 KT 1 1 10 1 10 2 10 3 
Photon Energy [GeV] 

FIG. 3: The prediction of our (i) best-fit model (thick lower 
green) and (ii) lighter DM mass (thin upper green), for the 
diffuse photon spectrum from the Inner Galaxy. The data 
correspond to the preliminary FERMI measurement [38]. The 
background is approximated by a power-law fitted to the data 
between 3 and 16 GeV. We also plot the primary, ICS, and 
FSR contributions to the signal of our best-fit model. The 
prediction of our best-fit model is consistent with the data, 
although it is slightly lower in the last bins. This is due to 
slight tension discussed in the text, between the preliminary 
data and constraints from HESS measurements of Galactic 
Center photons and the electron cosmic ray spectrum. 



we fit our decaying model to the PAMELA, FERMI 
and HESS leptonic results, and present predictions 
for the photon and neutrino signals, [40 fit only to 
the diffuse gamma excess. The disagreement with 
regards to the ICS contribution may arise, in part, 
from different DM masses and different propagation 
models. It would be interesting to see if annihilat- 
ing models can provide consistent fits to both the 
leptonic and diffuse gamma excesses. 
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